Nanotechnology, Science and Applications Dovepress 

open access to scientific and medical research 

9 Open Access FullText Article ORIGINAL RESEARCH 

Hysteresis modeling in ballistic carbon 
nanotube field-effect transistors 



This article was published in the following Dove Press journal: 
Nanotechnology, Science and Applications 
8 July 2014 

Number of times this article has been viewed 



Yian Liu' 

Mateus S Moura^ 

Ademir J Costa^'^ 

Luiz Alberto L de Almeida'' 

Makarand Paranjape' 

Marcio Fontana^ 

'Department of Physics, Georgetown 
University, Washington, DC, USA; 
^Department of Electrical Engineering 
Federal University of Bahia, Salvador, 
Brazil; 'Federal Institute of Bahia, 
Santo Amaro, Brazil; ""Engineering, 
Modeling and Applied Social Sciences 
Center, Federal University of ABC, 
Santo Andre, Brazil 



Correspondence: Marcio Fontana 
Department of Electrical Engineering, 
Federal University of Bahia, Salvador, 
40210-630, Brazil 
Email mfontana@ufba.br 



Abstract: Theoretical models are adapted to describe the hysteresis effects seen in the electrical 
characteristics of carbon nanotube field-effect transistors. The ballistic transport model describes 
the contributions of conduction energy sub-bands over carbon nanotube field-effect transistor 
drain current as a function of drain-source and gate-source voltages as well as other physical 
parameters of the device. The limiting-loop proximity model, originally developed to under- 
stand magnetic hysteresis, is also utilized in this work. The curves obtained from our developed 
model corroborate well with the experimentally derived hysteretic behavior of the transistors. 
Modeling the hysteresis behavior will enable designers to reliably use these effects in both 
analog and memory applications. 

Keywords: ballistic transport, nanoscale device, solid-state device 

Introduction 

Recent advances in nanotechnology have enabled the development and use of carbon 
nanotubes as field-effect transistors (CNTFETs). These devices have advantages over 
traditional metal oxide semiconductor transistors, including higher circuit integration,' 
increased gain, and higher mobility.^ The device consists of a single-wall (SWCNT) or 
multi-wall (MWCNT) carbon nanotube acting as the channel of a field-effect transistor, 
connected to two metal source and drain electrodes, and has its drain current controlled 
by a gate electrode. In general, a positive gate voltage shifts the Fermi level of the 
carbon nanotube toward the conduction band while a negative gate voltage shifts the 
Fermi level toward the valence band.' 

Hysteresis is a peculiar effect in CNTFETs and is observed when the gate voltage 
is swept in the forward and reverse directions, leading to a hysteretic effect that is 
difficult to predictably analyze due to the different nonlinear effects that are present. 
Hysteresis is typically related to charge injection from the carbon nanotube into the 
gate oxide dielectric, particularly at large gate bias, where the charges then become 
confined until the gate polarity reverses as a part of the voltage sweep.* Mathemati- 
cal models have been developed to characterize hysteresis behavior, such as the 
limiting-loop proximity (L^P),^ Preisach,'' and Jiles-Atherton' models. However, 
using only these models individually has not generated accurate approximations 
for CNTFET drain current versus gate voltage or /j^-K^ electrical characteristics 
primarily because they do not take into account the specific properties of the nano- 
tube. One theoretical model often used to describe CNTFETs was proposed,* based 
on the physical effects of ballistic conduction, wherein electrons in the nanotube are 
not affected by scattering.' Ballistic transport has nearly been achieved in CNTFET 
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modeling by the use of proper contact materials such as 
palladium and gold, thereby reducing transport barriers 
and approaching the theoretical conductance limit of the 
device (Go^»4qVh)."' 

The main objective of this work is to develop a reliable 
mathematical model by which to determine drain current by 
considering the physical basis of CNTFET conduction, which 
includes the hysteresis effect. Positive results were obtained, 
considering that some fitting parameters required estimation 
while instrument accuracy affected some of the obtained 
values. The results demonstrate that hysteretic effects can 
indeed be modeled and can potentially be used in a controlled 
manner for a variety of applications. 

In this paper we describe the ballistic theory pertain- 
ing to carbon nanotube transistor behavior, including the 
relationship between V^, V^^, /^^, and the transfer char- 
acteristics of the CNTFET. Adaptations to the L^P model 
to account for hysteresis, based on ballistic theory, are 
presented. We will also show the experimental procedures, 
methods, and materials used in order to obtain the device 
and data; explain the methodology used for the acquisition 
of the model parameters; present the results for different 
experiments and compare the data with simulations from 
the adapted model; and express conclusions about these 
comparisons. 

Model description 

Due to its dimensions, a carbon nanotube can be considered a 
quasi one-dimensional element, whose diameter and chirality 
directly affect its properties. There is no boundary scattering 
in pristine nanotubes due to the lack of surface bonds and a 
cylindrical structure. Thus, for long-length nanotubes, quasi- 
ballistic transport can be observed." '^ 

For a CNTFET, when contact is made between the carbon 
nanotube and the metal electrodes, a potential, or Schottky 
barrier is formed at the metal-carbon nanotube junction. 
This barrier has a height related to the work function of the 
metal electrode. Since gold was used as the contact mate- 
rial in this study, having a large work function, the Schottky 
barrier height became small. Hence, in this paper, we only 
consider the ballistic transport within the carbon nanotube 
channel and assume that the gate voltage only modulates the 
energy bands of the carbon nanotube channel between the 
source and drain.* 

The channel potential can be described as: 



(3, =Ke, for Fo <A1 
^,=Ko-t^(Fo-Al),forKo>Al 



where cp^ is the channel potential, is the gate bias, a is the 
slope of the drain current curve, and Al is the equilibrium 
sub-band minima of the first sub-band. 

In order to remove the discontinuity at Al , the following 
equation is proposed:'' 



a{V^-Al) + ^{aiV^-Al)f + e' 



(2) 



with £ serving as a smoothing parameter (typically around 
0.5 V). 

The slope of the characteristic curve, a, is expressed as a 
polynomial of V^^ such that a= 0C(^+ oc^ V^^ + V^^, where 
a^, ffj, and are dependent on both CNTFET diameter, d, 
and gate oxide thickness, t^^. The drain current, 7^^, can now 
be determined by the simplified expression below: 



/,,=^[ln(l + .^^)-ln(l 



+ e 



(3) 



where 



(4) 



and 



(5) 



thermal voltage 



q 



(6) 



(1) 



where V^^ is the threshold voltage, k is the Boltzmann con- 
stant, q the electron charge, and 7 the temperature. 

Hysteretic effects cause a shift of the threshold volt- 
age, which may be considered as the first sub-band, Aj. 
Experimentally, a sweep to large positive values of 
creates a positive shift in the threshold voltage,'^"" and 

does not affect the shift or the position of the curves. 
According to the variables set forth in the L^P model,' hys- 
teresis is characterized by a horizontal shift between the two 
curves, denoted as w, which designates the magnitude of 
shift around a central voltage, V^, along the horizontal axis. 
An auxiliary variable, 5, can be used to identify hysteresis 
orientation. These variables are experimentally obtained. 
In the case of CNTFETs, without these parameters shown 
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above, the current model cannot show hysteresis. Adapting 
the threshold voltage to L^P model variables, these effects 
are evidenced. 

Applying these definitions: 



(7) 



with 



+1, for — ^>0 
dt 

-1, for — ^ < 0 
dt 



(8) 



Although carbon nanotubes are described as having bal- 
listic behavior, the contact resistance between the electrode 
and the carbon nanotube'^ affects the on-state current level 
because the drain-source voltage in the nanotube is just a 
fraction of F^^ applied to the transistor: 



R V 

^^CN ' SD 



(9) 



where R^^ is the intrinsic resistance for the carbon nano- 
tube, which is determined as h/4q^ (approximately 6.4 kQ),^ 
while Rjj and R^ are the drain and source contact resistances, 
respectively. 

Figure 1 illustrates the conceptual procedure for CNT- 
FET modeling, where w and are obtained from the L^P 
model while the other variables arise from the ballistic 
model. 



Experimental procedure 

Silicon wafers (p-type, orientation <100>, 500-550 jim 
thickness, 0.001-0.005 Q cm resistivity) produced by 
NOVA Electronic Materials Ltd (Flower Mound, TX, USA) 
were used in this work. Thin silicon dioxide (SiO^) layers 
were grown on the wafer by a typical dry-wet-dry oxida- 
tion process at a temperature of 1 , 1 00°C for 1 0 minutes, 70 
minutes, and 10 minutes, respectively. Prior to the catalyst 
spinning procedure, the silicon wafer was cut into samples 
measuring 10x10 mm, which were ultrasonically degreased 
in trichloroethylene (C2HCI3), acetone ((CH3)2CO), and 
isopropyl alcohol (C^Hfi), rinsed in deionized water, and 
dried in nitrogen. A solution of ferric nitrate nonahydrate 
(1.6 mg) + precursor Mo02(acetylacetonate ligand)^ 
complex (0.5 mg) + aluminum oxide (15 mg) dissolved 
in methanol (20 mL) was used as the catalyst for carbon 
nanotube growth. Catalyst islands are patterned, using elec- 
tron beam lithography and lift-off on boron-doped silicon 
substrates with 300 nm SiO^ grown by thermal dry oxida- 
tion. SWCNTs are grown from these catalyst islands using 
chemical vapor deposition with a methane and hydrogen 
gas mixture (total flow 60 seem) in a 750°C-900°C-750°C 
three-zone ftirnace set. The SWCNT samples were grown 
with methane flow (32 seem) and hydrogen flow (28 seem). 
After imaging as-grown nanotubes with a field-emission 
scanning electron microscope, source/drain electrodes are 
patterned by electron beam lithography. The width of the 
source and drain metal electrodes are 1 |J,m with a 1.5 jim 
separation. The metals chromium and gold (thicknesses of 
30 A/1 ,000 A) are deposited by sputtering and followed by 
lift-off in acetone. 
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Figure I Flowchart for carbon nanotube field-effect transistor modeling. 
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Surface morphology of the grown layers was examined 
by field-emission scanning electron microscopy using a Zeiss 
SUPRA 55-VP FESEM microscope (Carl Zeiss Meditec AG, 
Jena, Germany). Micro-Raman spectroscopy was carried out 
at room temperature using a Renishaw Raman microscope 
(inVia; Renishaw, Wotton-under-Edge, UK) employing 
the output of an Ar+ laser (20 mW power) for excitation 
at X^514.5 nm. A Lab VIEW 2013 SPl software program 
(National Instruments Corporation, Austin, TX, USA) was 
implemented in an IEEE-488 environment, using a computer 
to control the Model 4 1406 pA meter/DC voltage source 
(Hewlett-Packard Inc., Palo Alto, CA, USA), and connected 
to a four-probe micromanipulator system for characterization 
of device parameters. All devices were measured at room 
temperature (300 K) in a controlled pressure (1 atm + 0.12 
kPa in H^O) and humidity (30%) environment. 

Model parameters identification 

For the determination of the L^P model parameters w and 
V^, we used graphic analysis based on the simulation results. 
Figure 2 describes the link between these parameters and 
their visual representation in the theoretical model, w is 
the horizontal distance between the two curves, obtained 
by tracing two lines in parallel with the inclination of the 
curves in their transition regions. Thereafter, is obtained 
by determining the horizontal midpoint of w at the vertical 
point located in the off-state saturation current value. 

Although a^, a^, and are inherent to the physical 
characteristics of the CNTFET, as mentioned in the Model 
Description section, there is no known mathematical relation 
that could allow the prior acquisition of these parameters. 
Thus, we considered each transistor as a stand-alone device. 



Saturation 




-i/_(max) 



+i/g(max) 



IVlidpoint: V 



Figure 2 /jd-V^- curves with hysteresis and obtainment of parameters w and V^. 
Notes: /^^ is the drain current: ^ ^ is the gate voltage; w represents the horizontal 
shift betv^een the tv^o curves; and represents the horizontal midpoint of w at the 
off-state saturation. 



obtaining the slope of the curve in the transition region for 
each device and drain-source voltage by the minimization 
of deviations between the data and the model curves in this 
region. In this way, with at least three measurements of a 
for different values of V^^, it is possible to determine a^, a^, 
and for a device by a linear system. 

Results and discussion 

Figure 3A shows a long and straight single semiconducting 
nanotube contacted by two gold electrodes to form a CNT- 
FET, and Figure 3B illustrates the typical hysteresis behavior 
of the nanotube transistor for drain current versus gate voltage 
(with Kj^=0.02 V). Data were taken at room temperature and 
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Figure 3 Fabrication and characterization of CNTFET. 

Notes: (A) Typical field-emission scanning electron microscope morphology 
of the carbon nanotube field-effect transistor with 300 nm silicon dioxide grown by 
thermal dry oxidation and source and drain metal electrodes (chromium and gold 
with thicknesses of 30 A/ 1, 000 A). (B) Typical hysteresis behavior of the carbon 
nanotube field-effect transistor at room temperature. Variables used are the drain 
current {^^^ and the gate voltage (V^.). 

Abbreviation: CNTFET, carbon nanotubes as field-effect transistors. 
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the device is operated with gate voltages between -10 V< 
<+10 y in the asymmetric hysteretic transition region. 
This result corroborates well with the data from a previous 
work," in which the CNTFET also showed that the device 
conductance was increased for negative gate voltages. We 
used the ballistic transport model' and L^P modeP for hys- 
teresis modeling of CNTFETs. Due to the effect complexity, 
this modeling is susceptible to nonlinearities, so we have 
incorporated effects in our model to include energy bands, 
memory phenomena, minor loop accommodation, and 
stabilization. 

Figure 4A and B, respectively, show typical Raman spec- 
trum and typical radial breathing mode peaks ranging from 
100 to 300 cm"' that were used to estimate the diameter of the 
SWCNT^" We observed two typical SWCNT peaks located at 
1,350 cm-' (D-band) and 1,590 cm"' (G-band). Table 1 shows 
details of the diameter distribution of S WCNTs synthesized in 
the furnace. The diameter distribution of the carbon nanotubes 
ranged between 1 .0 nm and 2.2 nm depending on the different 
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Figure 4 Raman spectra of the SWCNT. 

Notes: (A) Typical Raman spectrum of the single-wall carbon nanotube. (B) Typical 
radial breathing mode peaks which were used to estimate the diameter of the single- 
wall carbon nanotube. 

Abbreviation: SWCNT, single-wall carbon nanotube. 



Table 1 The diameter distribution of single 


-wall carbon nanotubes 


synthesized 


in the furnace 




KDn Da.ncis 




Relative intensity 


(cm ) 


(nm) 


Of the KDrl 


124 


2.2 


Weak 


133 


2.1 


Medium 


148 


1.8 


Weak 


167 


1.6 


Medium 


189 


1.4 


Medium 


198 


1.3 


Weak 


208 


1.2 


Weak 


259 


1.0 


Strong 



Abbreviation: RBM, Radial Breathing Mode. 



argon concentrations. In order to determine the carbon nanotube 
diameter, the Raman shift intensities were utilized and, based 
on the values in the table, we determined the diameter to be 
1.0 nm, so that Al was calculated as Al~0.45/d (nm) eV^' 
Contact resistance was measured in saturation current, and we 
obtained approximately 100 kQ. in this device. 

Simulation was performed at room temperature, using 
MATLAB (Math Works, Natick, MA, USA), on the device 
for different values of V^^ and different variances. For the 
device, ^„^=300 nm and we used fitting parameters aj=0.992, 
a =-0.05 V-' and a=0.024V-\ 

Figure 5A-C display curve modeling for between -2.5 
and +2.5 Y, -5.0 and +5.0 V, and -10.0 and +10.0 V, respec- 
tively, and for V^^ equal to 0.03 V, 0.02 V, and 0.01 V The 
results highlight how the gate-source voltage sweep affects 
hysteresis parameters and w. A larger sweep enlarges w 
and shifts toward negative values of V^, while changing 
does not affect these parameters. It is noted that the bottom 
curves of the model do not reach the same value of the upper 
curves on negative voltage limit, and the reason we suggest 
is that the devices did not reach their maximum saturation 
current, and, for this, the negative voltage limit should be 
lower. Despite this observation and some inaccuracies due to 
measurement and parameter deviation, a good fitting could 
be observed in all graphs. 

This work makes some contributions regarding CNT- 
FET modeling: first, simulations with higher biases of 
were performed, while several articles modeling CNTFET 
devices used a smaller gate voltage sweep (eg, 0 to +0.05 
volts) in simulations, in comparison to our work, and second, 
the description of hysteresis effect by using L^P parameters 
showed that hysteretic loop width and shift is sensitive to 
gate voltage bias. 

In order to evaluate the accuracy of the proposed model, 
three other devices were submitted to measures, for between 
-10 V and +10 V, with the same procedures and conditions 
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Table 2 Comparison of parameters and relative error among 
tested devices 




B 300 

250' 
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Figure 5 Comparison between experimental (red dot) and model (line) curves, for 
Vj„=0.03 V, 0.02 V, and 0.01 V. 

Notes: ranged from (A) -10 V to +10 V; (B) -5 V to +5 V; and (C) -2.5 V to 
+2.5 V. The variables used are the drain current {\^^, the gate voltage (VJ and the 
source/drain voltage (V^^). 

described earlier in the "Experimental procedure" section. 
Subsequently, the experimental results were employed in 
simulations for the modeling of the devices. Table 2 describes 
the parameters used for the devices and the relative error of 
the comparison between the experimental and model curves. 
It is noted that, with the results shown in Table 2, the model 
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«0 




a. 
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^c 


Relative 


number 


(V) 




(V-') 


(V-^) 






error (%) 
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0.01 


0.988 


-0.03 


2.00 


6.8 


3.4 


4.8 




0.02 
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0.03 












4.8 
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0.01 


0.990 


-0.01 


I.IO 


7.0 


2.0 


4.8 




0.02 












5.5 




0.03 












5.8 
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0.01 


0.990 


-0.04 


0.15 


6.4 


3.4 


5.7 




0.02 












4.8 




0.03 












2.4 



Notes: V^^ is the voltage between the drain and the source of the CNTFET; a^, a,, 
and are the fitting parameters specified in the model description: w and are 
parameters referent to the hysteresis effect in the device, explained in the model 
description as well. The relative error is between the modeled curve and the curve 
that represents the real data acquired, for each case. 



is applicable to other devices with different characteristics 
and behavior. In spite of the expected relative errors, since 
the measurement accuracy and model approximations might 
interfere in the modeling, these differences do not attenuate 
the effectiveness of the proposed work. 

Conclusion 

In this work, we successfully modeled 1sd~^g curves 
including hysteresis effects, for different voltage values 
and parameters, as a result of mathematical adaptations 
over the main expression, respecting the physical concepts 
of ballistic theory. The L^P model contributed to the char- 
acterization, as we could use its parameters to distinguish 
the two curves caused by hysteresis. The modeling brings 
detailed information about the device behavior, based on 
its characteristics. 

Currently, we are studying improvements in the char- 
acterization process. One improvement that we intend to 
demonstrate in a future work is the prior identification of 
the hysteretic parameters without using graphic methods. 
Nonetheless, the significant contribution presented in this 
article is that, given a device with its parameters, it is possible 
to mathematically localize a bias point over the experimental 
curves within a desired region of operation, thereby permit- 
ting accurate control of the CNTFET for specific applications, 
such as nonvolatile memory devices. 
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